Anomalous mass enhancement in strongly-correlated quantum wells 
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Using dynamical-mean-field theory, we investigate the electronic properties of quantum wells 
consisting of a t^-electron system with strong correlations. The special focus is on the subband 
structure of such quantum wells. The effective mass is found to increase with increase in the value 
of the bottom of the subband, i.e., decrease in the subband occupation number. This is due to 
the combination of Coulomb repulsion, whose effect is enhanced on surface layers, and longer-range 
hoppings. We discuss the implication of these results for the recent angle-resolved photoemission 
experiment on SrV03 thin films. 

PACS numbers: 73.21.-b, 71.10.-w 



Two-dimensional electron gases (2DEGs) realized in a 
variety of oxide interfaces have been attracting significant 
interests*^ In particular, electronics utilizing oxides with 
strong correlations would benefit from their rich phase 
behaviors.— For example, control of the band structure of 
2DEGs in transition-metal oxides has been proposed as 
a way to create non-cupratc high-T c superconductivity^ 
Yet, the realization of metallic behavior in few-unit-cell 
thick oxides remains challenging^^ 

Two-dimensional metallic behavior in confined geom- 
etry, i.e., in quantum wells (QWs), has been stud- 
ied for conventional metals. Reconstructed band struc- 
tures or subband dispersion relations in QWs of Ag 
thin films have been confirmed using photoemission 
spectroscopy^ The subband dispersion of 2DEGs re- 
alized on the surface of a band insulator SrTi03 was 
also observed using angle-resolved photoemission spec- 
troscopy (ARPES)^^ More recently, Yoshimatsu and 
coworkers have performed ARPES measurements on 
QWs in thin films of the correlated metal SrVOs^i The 
subband structures realized in such QWs can be ex- 
plained reasonably well using a simple tight-binding-type 
description. However, the effective mass of such sub- 
bands was found to increase with decreasing binding en- 
ergy of the subband. Since this trend is opposite to what 
we expect based on the bulk behavior, i.e., the effective 
mass is reduced with decreasing binding energy and de- 
creasing band occupancy, the origin of such an anomalous 
mass enhancement remains to be understood. 

In this Rapid Communications, we analyze model 
QWs consisting of a t\ g electron system as experimen- 
tally considered by Yoshimatsu et al. We employ layer 
dynamical-mean- field theory (DMFT) with the exact- 
diagonalization impurity solver] 12 ' 13 In correlated QWs, 
a smaller coordination number on surface layers induces 
larger mass enhancement than in the bulk region^ - — 
This brings about the anomalous subband-depcndent 
mass enhancement; the effective mass becomes larger 
with decreasing subband binding energy or depopulation 
of the subband. With the additional effect coming from 
the longer-range hopping, the subband-dependent mass 
enhancement was found to increase dramatically. We ar- 
gue that the anomalous mass enhancement reported for 



thin films of SrVC>3 is due to strong correlations and long- 
range transfer integral. 

We consider the three-band Hubbard model involving 
t2 g electrons, H = .ffband + Sr^ioc( r )- The first term 
describes the non-interacting part of the system as 

-ffband = - ^2 ^2 trr'dr Ta dr'Tcr, (1) 

T,a YX' 

where dr T a stands for the annihilation operator for an 
electron at site r in orbital r with spin a, and i£ r , is 
the hopping integral between orbital r at sites r and r'. 
For the orbital r = xy, we take the nearest neighbor 
hoppings i£ r , = t n for r' = r ± x(y) and t£ r , = t$ for 
r' = r±z, and the second-neighbor hopping t£ r , = t a i for 
r' = r ± x ± y. Here, x(y, z) is the unit vector along the 
x(y, z) direction. The hopping parameters for the other 
orbitals are given by interchanging the coordinates x, y, 
and z accordingly. Parameter values are taken from den- 
sity functional theory results as t v = 0.281, ts = 0.033, 
and t a > — 0.096 (all in eV)^ Hi oc describes the local 
interaction as 

t,t',t n T ■ cr 

r // \cr,a' 

(2) 

Here, the site index r is suppressed for simplicity, and \i 
is the chemical potential. Since we consider ti g electron 
systems, we assume the relation U = U' + 2J, where 
U = U TTTT (intraorbital Coulomb), U 1 = U TT ' TT ' (in- 
terorbital Coulomb), J = U TT T T (interorbital exchange) 
_ jjttt t (interorbital pair transfer) for r ^ t', and 
other components are absent As in Ref. [Tl|, we con- 
sider QWs in which a finite number of correlated layers 
stack along the z direction with the open-boundary con- 
dition and the periodic-boundary condition along the x 
and y directions. 

Before going into the detailed analysis taking into 
account the correlation effects, let us first discuss the 
low-energy electronic behavior focusing on quasiparticle 
bands. For this purpose, we consider the following effcc- 
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FIG. 1: (Color online) (a) Dispersion relation as a function of 
momentum for a non-interacting 5-ML-thick quantum well. 
Crosses indicate E\f a (binding energy at k = times —1). 
(b) Three-dimensional plot of subband quasiparticle weight 
Zg z , El za , and QW thickness. For each QW, Z»* is defined 
at E^_°^ < 0. When projected on the left (right) vertical 
plane, Z^f is given as a function of Ey^J* (QW thickness). 



tive one-dimensional Schrodinger cquatio: 



,21 



z,z'±l 



= ^kXk> ( 3 ) 

where EJ(w) is the electron self-energy at orbital r in 
layer z computed in layer DMFT, and is the layer- 
dependent quasiparticle weight defined by Z T Z = {1 — 
Re9 CJ E^(o;)| LJ= o}~ 1 - ^j. is the in-plane dispersion and 
is the out-of-planc hopping clement for orbital r with in- 
plane momentum k = (k x ,k y ). For t = yz, these are 
explicitly given by = — 2t ff cos k y — 2t a i cosk x and 
fQ = —tir — 2t§ cos k y . a labels the subband with the 
eigenfunction <p r ^ in increasing order of the subband en- 
ergy . As an example, the energy eigenvalue 
for 5-ML-thick non- interacting [Z T Z = 1) QW is plotted 
in Fig. Q] (a) (ML indicates monolayer). We notice that 
subbands originating from yz{xz) orbitals are not paral- 
lel, while xy subbands are. This is because the second- 
neighbor hopping between neighboring layers induces k 
dependence in the out-of-plane hopping i^ xz \ For or 



ijr — 2ts cos k y ( x ) and, therefore, the 



bital yz(xz), 

subband separation becomes large when k y ( x \ approaches 
0. As a result, the Fermi velocity of high-energy (less- 
populated) bands becomes small as if the effective mass 
is enhanced. 



The low-energy electronic behaviors of correlated QWs 
are governed by the quasiparticle subbands. The cor- 
relation effects enter as the quasiparticle weight of the 
subband. Using the solution of Eq. ((3]), the subband- 
dependent quasiparticle weight is given b y 21 ' 22 



^=E^ r kk=k" 



(4) 



From Eq. it is clear that the subband quasiparti- 
cle weight becomes unity in the absence of correlations, 
i.e., Z T Z — 1 leads to Szl^k-k" I = ^ (normalization 
of the quasiparticle eigenfunction). Another important 
quantity is the effective quasiparticle weight defined by 



ko 



(5) 



k=kl 



Here, k™ is the Fermi momentum for the crth subband, 
and e 



_ ~ y z 

£ kk, - £ k 



k is the bulk dispersion. For t — yz, we have 

fQ — 2i 7r cosfc z . Thus, Z T a measures the 
change in the Fermi velocity with respect to its bulk 
value. In Ref. [n], Z T a was used to discuss the mass en- 
hancement. 

Because of the k dependence of the interlayer hopping 
matrix tT , Z can be smaller than unity even without cor- 
relations. Figure [T] (b) summarizes the results for E^ a 

and 1/Z for non-interacting QWs with thickness varied 
from 4 to 10. As the QW becomes thin, increases 
and the number of occupied subbands is reduced (see the 
basal plane). For a 5-ML-thick QW, 1/Z^ 2 is projected 
in the left vertical plane thus is shown as a function of 
^k— o (binding energy times —1). As E^ a approaches 
0, l/Z^ z is increased. This trend can be seen in all QWs 
studied (see the projection of l/Z^f on the right verti- 
cal plane). A similar trend was reported experimentally. 
However, \jZ T a is enhanced from ~ 1 to ~ 1.7, so is at 
most 70 %. Therefore, the band effect alone does not ac- 
count for the large mass enhancement reported in Ref. ITlL 
where varies from ~ 1.7 to ~ 4.5. The experimen- 

tal enhancement in l/Z^ is nearly 300 %, and at 
the largest binding energy is already ~ 70 % larger than 
the band mass. These results indicate the influence of 
the correlation effects. 

In order to see the effect of correlations rather quan- 
titatively, here we employ layer DMFT, whose self- 
consistency condition is closed by^ ~ 16 ' 21 



(6) 



Here, G T Z is the local Green's function on layer z, and 
the lattice Green's function matrix on the right-hand 
side is given as a function of k and the z-axis coordi- 
nate as G(k,w) = [(lj + — Hband(k) — S(k,w)] 




Thickness (ML) 

FIG. 2: (Color online) (a) Orbitally-resolved spectral function 
as a function of momentum and frequency for the interacting 
5-ML-thick QW with U = 4 and J = 0.5 eV. Dotted lines 
indicate the solution of Eq. ([3]), , obtained using the layer 
DMFT result for the quasiparticle weight Z T Z . Note that yz 
and xz bands are symmetric with respect to X Y. (b) 
Subband quasiparticle weights Z^ z and Z^ as functions of 
E^. (c) E^ a , (d) quasiparticle weight Zg", and (e) effec- 
tive quasiparticle weight Z^ z as functions of the thickness of 
QWs. Plots (b)-(e) are generated using the method displayed 
in Fig. [T](b). See also Ref.[H (Fig. SI). 

The hopping matrix -ffband(k) is given by an in-plane 
Fourier transformation of -ffband as -ffband(k) = (s^5 ZyZ > + 
t^S z-z >±i)S T-T > . The self-energy matrix is approximated 

as E(k,w) = T<1(uj)5 z , z 'S t . t i . The local self-energy is ob- 
tained by solving the effective impurity model defined 
by the local interaction term coupled with an effective 
medium. In this study, we use the exact diagonalization 
impurity solver with the Arnoldi algorithmi 23 i 24 Here, the 
effective medium is approximated as a finite number of 
bath sites, and the impurity Hamiltonian is given by 

-ffimp = H loc + ^ £irc\ T(T CiT<r + ^ {^ir^ iTG d r( , + h.cA . 
i,T,a i.r.a 

(7) 

Ci Ta is the annihilation operator of an electron at the zth 
bath site with potential £ JT and hybridization strength 
with the impurity orbital r denoted by Vi T . Because of 
the exponentially growing Hilbert space with respect to 
the numbers of orbitals and electrons, we consider two 
bath sites per correlated orbital, i.e., i = 1,2. In our 
numerical simulations, we use temperature T = 10~ 2 eV 
to retain low-energy states with Boltzmann factors larger 
than 10~ 6 and consider only paramagnetic solutions. 

Figure [5] (a) shows the results for the orbitally-resolved 
spectral function A T (k 7 ui) = — i J2 Z ImGJ z (k,a;) as well 
as ET a as dotted lines for a 5-ML-thick interacting QW 
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FIG. 3: (Color online) (a) Orbitally-resolved occupation num- 
ber as a function of position z for a non-interacting 10-ML- 
thick QW. (b) Same as (a) for the interacting model and (c) 
the local quasiparticle weight as a function of z. Layer DMFT 
with (7 = 4 and J = 0.5 eV is used for (b) and (c). Inset: 
Quasiparticle wave functions f v -\"_ a ■ 



with U = 4 and J = 0.5 (both in eV). For A T (k,uj), the 
self-energy is extrapolated to the real axis using the Pade 
approximation^ In comparison with the non-interacting 
case, the bandwidth is reduced by about 50 % due to 
correlations. We notice that E^ 1 reproduces the peak 
positions of A T {k 1 ui) fairly well. There are five sub-bands 
for both yz and xy, but those in the latter are indis- 
tinguishable because all subbands are located within the 
range of 2tg ~ 0.07 eV. Thus, we focus on yz subbands in 
the following discussion. Using the same procedure as in 
Fig.[T](b), we analyze E^ a ^ and the mass enhancements 

l/ZJ? and \jZl z . Figure El (b)shows plots of l/Z% z and 
l/Zv* as functions of E^" Q for a 5-ML-thick QW. At 

the largest binding energy, both 1/Z^ Z and l/Z^ 2 are 
about 2, the mass enhancement expected in the bulk re- 
gion. Figures [2] (c), [2] (d), and [2] (e) summarize E? z a , 

1/ZV Z , and l/Z^, respectively, for interacting QWs with 
thickness varied from 4 to 10. Although 1/Z^ Z shows an 
increase with increasing it is only from ~ 2 to 

~ 2.2. On the other hand, 1/Z„ Z shows a rather steep 
increase from ~ 2 to ~ 3.8, as reported experimentally^ 

Aside from the quantitative difference, non-interacting 
QWs and interacting QWs behave quite similarly. As a 
small but clear difference, some of interacting QWs have 
a larger number of occupied yz subbands, 7 and 9-ML- 
thick QWs. This is caused by the different orbital polar- 
ization. As shown in Figs.[3](a) and[3](b), non-interacting 
QWs have larger orbital polarization with smaller occu- 
pancy in the yz and xz orbitals on surface layers than in- 
teracting QWs. This is because these bands have a quasi- 
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FIG. 4: (Color online) Subband mass enhancement, (a) Effec- 
tive mass enhancement 1/Z^ Z for the non-interacting model 
as a function of E^_^. (b) Mass enhancement 1/Z% Z and (c) 

1/Z^ Z for the interacting model with [7 = 4 and J = 0.5 (eV). 
(d) 1 / Z^ z for the interacting model with the parameter values 
indicated. Gray bold lines are guides to the eye. 



one-dimensional character on the surface layers, with a 
reduced effective bandwidth. The average charge density 
also shows Fricdcl-type oscillatory behavior with respect 
to z. On the other hand, in the interacting case, the or- 
bital polarization and the charge redistribution are sig- 
nificantly suppressed because the charge susceptibility is 
suppressed near integer fillings. This behavior was found 
to be insensitive to the choice of the interaction strength 
as shown in Ref. [2^ (Fig. S2). Therefore, the effect of 
charge relaxation is expected to be small, in contrast to 
LaTiOs/SrTiOa heterostructures. 

Figure [3] (c) shows the position-dependent quasiparti- 
cle weight ZJ, and its inset the quasiparticle eigenfunc- 
tions for a yz electron at the Fermi level. Strong mass 
renormalization takes place in surface layers where the 
coordination number is smaller In the current case, 
ZZ, is <~ 0.43(0.51) on surface layers (in the bulk region at 
2 = 5); thus there is about 15 % stronger mass renormal- 
ization on the surface. This small difference comes from 
the fact that SrVC>3 is not so strongly correlated. With 
increasing U, Z T Z on surfaces are more strongly renor- 
malized [see Ref. HH, Fig. S3 (a)]. Since an eigenfunction 
with larger a has larger weight on the surface layers, 
the effective mass of such a subband is more strongly 
renormalizcd. But, the renormalization of Z^ z , up to 
~ 10 %, is smaller than that of Z T Z because of the in- 
terlaycr hybridization. The additional enhancement in 



1/Z^ Z is caused by the momentum dependent inter-layer 
hoppings, as discussed earlier. 

Figure S] summarizes the mass enhancement as a func- 
tion of the positioii_of the bottom of the subband. The 
enhancement in 1/ 1 Z"^ z is rather small for non-interacting 
QWs because it comes from the small hopping parame- 
ter t a i . The mass enhancement 1 jZ^ originating purely 
from the correlation effects also shows rathcr^small de- 
pendence on On the other hand, 1/Z^ Z shows 
strong -Ej^lg dependence because both the band and 
correlation effects arc included. The effective mass en- 
hancement 1/ZV Z somewhat depends on the correlation 
strength, as shown in Fig. @] (d). Comparison of l/Z'^ 2 
for different interaction strengths is presented in Ref. |25| 
[Figs. S3 (b)-S3 (c)]. With a reasonable parameter set, 
the experimentally reported mass enhancement can be 
scmiquantitativcly reproduced. 

We notice that the number of subbands is overesti- 
mated by ~ 1 for interacting QWs compared with the 
experimental observation^ A possible explanation for 
this discrepancy is that, in the experiment of Ref. [Ill, 
the surface layer is made of VO2, so that the symmetry 
and the valence state of surface V ions greatly deviate 
from those in the bulk. Also, we cannot exclude the pos- 
sibility of surface lattice relaxation by which conduction 
electrons are strongly localized on the surface layer. In 
these cases, the surface V sites would not contribute to 
the ARPES spectrum near the Fermi level as do those 
in the bulk. Detailed study including these effects would 
be necessary to fully understand the nature of SrV03 
QWs including the dimensional crossover and the metal- 
insulator transition£iii Yet, the present study provides a 
reasonable account for the anomalous mass enhancement 
reported for SrV03 thin films. 

Summarizing, using dynamical-mean-ficld theory, we 
investigated the electronic properties of correlated quan- 
tum wells consisting of a i^-electron system. The special 
focus is on the subband structure of such quantum wells. 
The subband effective mass was found to increase with 
decreasing band occupancy as reported for SrV03 thin 
films. The present theory provides a reasonable account 
for this observation as the combined effect of Coulomb 
repulsion, whose effect is enhanced on surface layers, and 
longer-range hoppings. Inclusion of these two effects is 
essential to correctly interpret experimental observations. 
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Supplementary material: Anomalous 
mass enhancement in strongly-correlated 
quantum wells 

For interacting quantum wells (QWs), each subband 
has the quasiparticle weight Z T a defined in Eq. (4) and 



the effective quasiparticle weight Z T a defined in Eq. (5). 
Figure SI shows three-dimensional plots of (a) Z"^ and 




FIG. SI: Three-dimensional plot of Eg" < 0, QW thickness, 
and Z^ z (a) and Z% z (b). Here, the interaction strength is 
taken as U = 4 and J = 0.5 eV. For each QW, Z^ z is defined 
at Eg™ < 0. When projected on the left (right) vertical 

plane, Z^ z is given as a function of Eg_^ (QW thickness). 

with C/ = 4 and J = 0.5 cV. On the basal planes, Eg" 

is plotted. When Z^(Z^) is projected in the left (right) 
vertical plane, it is given as a function of E? z " (QW 

It— 

thickness). 

Figure S2 summarizes the subband effective mass 
l/Z^ z as a function of E? z " , i.e., the binding energy 

it— 

times — 1; [7 = 3 and J = 0.5 (top panel), U = 4 and 
J = 0.5 [middle panel, same as Fig. 4 (b)], U = 5 and 
J = 0.5 eV (bottom panel). With increasing U, 1/Z% Z 
is increased. At the same time, l/^ z at larger Eg" is 
enhanced more strongly. 

Figure S3 (a) shows the quasiparticle weight Z T Z as 
a function of layer index z for 10-ML-thick QWs with 
different interaction strengths. With the increase of U, 
Z T Z on surface layers is more strongly rcnormalized. 

Figures S3 (b)-(c) show the orbitally resolved charge 
density as a function of layer index z for 10-ML-thick 
interacting QWs. The charge redistribution, deviation 
of the average density from 1/3, is suppressed with the 
increase of U. 
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FIG. S2: Subband quasiparticle weight Z% z as a function of 
^k— o ^ e binding ener Sy times —1) with different interaction 
strengths indicated. 
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FIG. S3: (a) Quasiparticle weight ZJ with r = xy,xz,yz as 
a function of layer index z for 10-ML-thick QWs with differ- 
ent interaction strengths indicated. With increasing U, ZZ is 
renormalized, and the renormalization becomes stronger on 
surface layers at z = 1 and 10. Orbitally resolved charge den- 
sity as a function of z for 10-ML-thick QWs with (7 = 3 and 
J = 0.5 eV (b), (7 = 4 and J = 0.5 eV (c), and U = 5 and 
J = 0.5 eV (d). The average charge density is also shown. 
Gray dotted lines indicate 1/3. 



